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Abstract 

This paper is concerned with an inverse obstacle problem which employs the dy- 
namical scattering data of acoustic wave over a finite time interval. The unknown 
obstacle is assumed to be sound-soft one. The governing equation of the wave is 
given by the classical wave equation. The wave is generated by the initial data lo- 
calized outside the obstacle and observed over a finite time interval at a place which 
is not necessary the same as the support of the initial data. The observed data 
are the so-called bistatic data. In this paper, an enclosure method which employs 
the bistatic data and is based on two main analytical formulae, is developed. The 
first one enables us to extract the maximum spheroid with focal points at the center 
of the support of the initial data and that of the observation points whose exte- 
rior encloses the unknown obstacle of general shape. The second one, under some 
technical assumption for the obstacle including convexity as an example, indicates 
the deviation of the geometry of the boundary of the obstacle and the maximum 
spheroid at the contact points. Several implications of those two formulae are also 
given. In particular, a constructive proof of a uniqueness of a spherical obstacle 
using the bistatic data is given. 
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1 Introduction 

In this paper, we consider an inverse obstacle scattering problem for a sound-soft obstacle 
with dynamical data over a finite time interval. The governing equation of the wave is the 
classical wave equation. The wave as the solution is generated by the initial data whose 
support is localized at the outside of the obstacle and observed over a finite time interval 
on a different position from the support of the initial data. The observed data are the 
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so-called bistatic data. This is a simple mathematical model of the data collection process 
using an acoustic wave/electromagnetic wave such as, bistatic active sonar, radar, etc. 
See, e.g., [4] for the bistatic active sonar. The aim of this paper is to develop an enclosure 
method which employs the bistatic data. 

Let us describe a mathematical formulation of the problem. Let D be a nonempty 
bounded open subset of with C^-boundary such that \ D is connected. Let < 
T < oo. Let / e L^(R^) satisfy supp/ f] D — Let u — Uf{x,t) denote the weak 
solution of the following initial boundary value problem for the classical wave equation: 

d^u - Aii = in (R^ \D)x ]0, T[, 
u{x,0) = inR^ \;D, 

(LI) 

dtu{x,0) = f{x) mR^\D, 
u^O ondDx]0, T[. 

Here u denotes the unit outward normal to D on dD. The boundary condition for u in 
(LI) means that D is a, sound-soft obstacle. In this paper, T is always fixed. Thus, for 
our purpose the weak solution over the bounded interval ]0, T[ is appropriate. Since the 
notion of the weak solution for the wave equation is well established, we do not repeat 
the description here. Instead see [5] for the notion and also its use in [9, 11] for inverse 
obstacle scattering problems with dynamical data over a finite time interval. 

In this paper, we consider the following problem. 
Inverse Problem. Let B and B' be two known open balls centered at p e R^ and 
p' e R^ with radius r] and rj', respectively such that B H D — and B fl D = 0. Let 
Xb denote the characteristic function of B and set / = xb- Assume that D is unknown. 
Extract information about the location and shape of D from the data Uf{x, t) given at all 
X eB' eiudte ]0, T[. 

As far as the author knows, there is no result to this problem for general configulation of 
B and B'. This is the problem raised in [11] as an open problem related to the enclosure 
method itself. In particular, the problem contains the case when B (1 B' = (/} which 
corresponds to the case when the emitter and receiver are placed on different positions at 
a finite distance from the obstacle. Strictly speaking, we should call the data in this case 
the bistatic data, however, we include also the case B f] B' ^ $. 

In this paper, we develop an enclosure method with bistatic data. In short, the 
enclosure method aims at extracting a domain that encloses an unknown discontinuity, 
such as cavities, cracks, inclusions or obstacles. The idea of the enclosure method goes 
back to [7], in which the original enclosure method was developed by considering an 
inverse boundary value problem governed by the Laplace equation. In [8], an idea for 
the application of the enclosure method to the dynamical data coming from the heat or 
wave equations has been introduced. Now we have many applications of this enclosure 
method to inverse boundary value problems governed by the heat equations in [14, 15, 10], 
visco-elastic system of equations [13] and inverse obstacle scattering problems governed 
by the wave equations in [9, 11, 12]. 

We establish two main analytical formulae. The first one enables us to extract the 
maximum spheroid with focal points at the center of the support of the initial data and 
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that of the observation points whose exterior encloses the unknown obstacle of general 
shape. The appearence of the exterior of a spheroid as an enclosing domain is new since 
previous enclosing domains arc a half plane/space, sphere or its exterior, or cone. The 
formula shows us an effect of the bistatic data on the obtained information. See Theorem 
1.1 below. The second one, under some technical assumption for the obstacle including 
convexity as an example, indicates the deviation of the geometry of the boundary of the 
obstacle and the maximum spheroid at the contact points. This is also new. See Theorem 
1.3 below. And also we present several implications of those two formulae. In particular, 
we give a constructive proof of a uniqueness of a spherical obstacle using the bistatic data. 

1.1 Extracting the first refiection distance and its implication 

In this paper, given an arbitrary h e L^(R^), we denote by Vh the unique weak solution 
V e H^{Y{?) of 

(A - t'^)v + h{x) = in R^. (1.2) 

Vh has the expression 

1 [ Q-i'\^-y\ 

Vh{x) = Vh{x,T) = — rh{y)dy. (1.3) 

47r >/r3 \x — y\ 

Define 

rT 

Wf{x) = Wf{x, t) = e~'^^Uf{x, t)dt, x eK^\D, r > 0. 

t/ 

w = Wf satisfies 

{A-T^)w + f{x) =e-^^Ff{x,T) mR^\D, 

(1.4) 

w — ondD, 

where 

Ff{x, t) = dtUf{x, T) + TUf{x, T), xeR^\D. 

Since Ff{x,T) is unknown, it seems that the existence of such term in (1.4) hides the 
information about an unknown obstacle. However, the use of the enclosure method pre- 
sented below docs not make it a problem at all and enables us to extract the information 
about the obstacle provided T is sufficiently large and fixed. 

Let xb' denote the characteristic function of B' and set g = xb'- 
The results of this paper are concerned with the asymptotic behaviour of the indicator 
function: 



(/^^g - Wfg)dx = ^ Vgdx - j^^ Wfdx 



T I > _ 

/R3\D 

For the description of the results we prepare some notation. 
Define 

0(x; y,y')^\y-x\ + \x-y'\, (x, y, y') e x x 

This is the length of the broken path connecting y to x and x to y' which plays the central 
role in this paper. 

In this paper we denote the convex hull of the set F C R^ by [F] . 
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Theorem 1.1. Let [B U b'] f] dD ^ $ and T satisfy 

^> S^R ,.«oX^;?/'^')- (1-5) 

x€dD, yeoB, y'^oB' 

Then, there exists a tq > such that, for all r > tq, 

J^^^-{fVg-Wfg)dx > 

and the formula 

lim —log/ (fVo — Wfq)dx=— min d>(x:y,y') (1.6) 

Z5 valid. 

Note that 

SIr ,.«r,'^(^;^'^') ^ min - (77 + 77')- (1-7) 

x&dD,y&dB,y'&dB' x&oD 

See Appendix for the proof of (1.7). The quantity min^gg£)0(a;;p, p') coincides with the 
shortest length of the broken paths connecting p to a point q on dD and q to p', that 
is, the first reflection distance between p and q by D. (1.6) gives an extraction formula 
of mm^^QD ((){x; p, p') from Uf{x,t) given at all x & B' and t e]0, Tf. Formula (1.6) gives 
the method of carrying out calculation processing of the waveform mathematically, and 
extracting the first reflection distance. 

Define [p,p'] = {sp+{l — s)p' | < s < 1}. This is the straight line segment connecting 
the centers of B and B' and coincides with [{p,]?'}]. Since both p and p' are in \ D, 
[p,p'] n D = if and only if \p,p'] n = 0. 

We know that 

• min^ea£,0(a;;p,p') > \p - p'\; 

• if \p,p'] n dD = 0, then min^^Qj:) (p{x; p, p') > \p — p'\. 
Given c > \p — p'\ define 

Ec{p,p') = {x e I 4>{x;p,p') = c}. 

This is a spheroid with focal points p and p'. It is a compact surface of class C°°. 

Since D is contained in the exterior of spheroid Ec{p,p') with c = min^goc (l){x;p,p'), 
Theorems 1.1 gives us the largest spheroid with focal points p and p' whose exterior 
contains D using dynamical bistatic data -Uj on B'x ]0, T[. The appearance of the spheroid 
in the enclosure method is new and this is a decisive difference from the previous enclosure 
method. 

Therefore we obtain the information that there exists a point belonging to dD on the 
spheroid Ec{p,p') with c = min^-ga/) (f){x;p,p') calculated by formula (1.6). Thus, the next 
problem is: identify all the points belonging to dD on the spheroid. In order to describe 
the problem precisely we introduce the following notion. 
Definition 1.1. Let p and p' satisfy \p,p'] fl dD — 0. Define 

Asd(p,p') = {qedD\ (l){q;p,p') = mm(l){x;p,p')}. 

xedD 
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We call this the first reflector between p and p'. The points in the first reflector are called 
the first reflection points between p and p'. Note that h.QD{p-iP') can be an infinite set. 
One has the expression 

KdD{p,p') = dDf\E,{p,p') 

with c = min^g^ij Thus the problem becomes: identify all the first reflection 

points. 

Let uj G S"^. We denote by 5(0;; p, p', c) the length of the straight line segment con- 
necting p' and the unique point on Ec{p,p') H {p' + sou \ s > 0}. We have 

, c^-\p-p'\^ 
s{uj;p,p,c) - 



2{c — (jj ■ {p — p')} 

Note that u ■ {p — p') < c since c > \p — p'\. It is easy to see that the map 

S"^ 3 u ^ p' + s{u}; p, p', c)u} e 

is one-to-one and the image coincides with Ec{p,p'). 

Let < ri' < inf(^g52 s(a;;p,p', c). B is contained in the set of all x such that 
4>{x]p,p') < c, that is, the domain enclosed by Ec{p,p'). 

The following theorem says that all the first refiection points between p and p' together 
with the tangent planes can be extracted from a single set of the bistatic data. This 
exceeds the previous enclosure method and suggests that the information which contained 
in the bistatic data is quite rich. 

Theorem 1.2. Assume that c — min-rgg^i (f){x;p,p') is known. Let [B U B'] O D — 0. Fix 
< s <r]'. IfT satisfies 

T > sup min (f){x;p,p' + sou) — {rj + r]' — s), (1.8) 

then, one can exiract all q e Aqd(p,p') together with Vq from Uf on B'x]0, T[ with 
f = Xb- 

Remcirk 1.1. Note that 

sup min (j){x;p,p' + suj)-{ri + ri' -s)= sup "?^5o .H^'^V^v')- 

The reason is the same as that of the validity of (1.7). Thus the constraint on T is 
reasonable. 

Theorem 1.1 is a direct consequence of the following two estimates: there exist Hj G R, 
Cj > with j = 1,2 and Tq > which are independent of r such that, for all r > Tq, 



grmin,gao,^eas,j/'eas' 'Pi^;y,y') J ^^_{fVg - Wfg)dx < CiT^' (1.9) 
and 

C2T^^ < e^'^^-^9'^'y^9'^'y'^9B'H^;y,y') J ^^_{fvg - Wfg)dx. (1.10) 
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The proof of (1.9) proceeds along the same hne as the back-scattering data case {B — 
B') and is given in Section 2. The point that should be emphasized in the proof of 
Theorem 1.1 is (1.10) which is proved in Subsection 3.1. 

When B' = B, using the same technique as done for the sound-hard obstacle case in 
[11], we can prove (1.10) without difficulty. The technique therein does not depend on the 
boundary condition, however, heavily depends on the condition B — B'. In this paper, 
we take another way. It is based on the combination of the maximum principle for the 
modified Helmholtz equation in the domain \ D and a reflection across dD. It heavily 
depends on the speciality of the homogeneous Dirichlet boundary condition on dD. The 
idea goes back to the arguments done in the proofs of Theorem 3.6 and Lemma 3.7 in 
Lax-Phillips [18]. Note that, therein, a relationship between the support function and the 
so-called scattering kernel for a general sound-soft obstacle has been established. They 
used the arguments to obtain an estimate for the analytic continuation of the Fourier 
transform of the scattering kernel and then applied the Paley-Weiner theorem. We refer 
the reader to [19, 22] for several other results using the scattering kernel. 

1.2 Leading term of the indicator function and its implication 

(1.9) and (1.10) suggest that the following integral as r — > oo 



may behave as some power of r multiplied by a positive constant. The constant may 
contain some information about the geometry of the boundary of the obstacle at the 
points on dD that attain T^irix£dD,yedB,y'edB' (p{x;y,y'), i.e., the first reflection points 
between p and p'. 

If g e Aq£,{p,p'), then q G Ec{p,p') with c = min^^Qi:) (f){x; p, p') and the two tangent 
planes at q of dD and Ec{p,p') coincide. We denote by Sq{dD) and Sq{Ec{p,p')) the shape 
operators (or the Weingarten maps) at q with respect to Uq. Those are symmetric linear 
operators on the common tangent space at q of dD and Ec{p,p'). It is easy to see that 
Sq{Ec{p,p')) — Sq{dD) > as the quadratic form on the same tangent space at q (see 




(4.21)). 

Let Zf e H\R^\D) solve 



(A - T^)Zf = Ff{x,T) inR^ \ D, 



(1.11) 



Zf = ondD. 

It follows from (1.2) ior h — f and (1.4) that Wf has the form 

= f^/ + 4 + 



(1.12) 



where satisfies 



(A-T2)e5 = 0inR3\:D, 




Note that: since supp fr\D — $,Vf is smooth in a neighbourhood of D and thus, by elliptic 
regularity, we see that e° is smooth for x E 'R^ \ D. Moreover, note that e°(x) — > 
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as \x\ — > oo rapidly and uniformly with respect to x/|a;|. This is a combination of 
the uniqueness of the weak solution of (1-13) and a potential theoretic construction of 
the solution, see, e.g., [3, 20] for the approach and [15] for an application to an inverse 
problem for the heat equation. 

Given a; e define dgoix) = mfy^QD \y — x\. It is well known that there exists a 
positive constant So such that: given x e D/x E \ D with ddoix) < 25o there exists 
a unique q — q{x) be the boundary point on dD such that x — dgu{x)i'q ([6]). One 
may assume that both d^jji^x) and q{x) is for x G with dgn^x) < 2So; x e \ D 
with dQu{x) < 26o (Lemma 1 of Appendix in [6]). Note that i/g is the unit outer normal 
to dD at q. For x with dgD{x) < 25o define x^ = 2q{x) — x. 

Before describing our third result, we introduce a restriction on a class of obstacles 
which is satisfied with all convex obstacles. 

Definition 1.2. We say that D is admissible, if there exist positive constants C, S'{< 2So) 
and To such that, for all y & D with dgoiy) < S' and r > Tq 



I ej (?/'■) I <C f e-^l^-^lrfx. 

J B 



The following theorem gives an answer to the question raised above. 

Theorem 1.3. Let B and B' satisfy [S U n !D = 0. Let f ^ xb and g ^ xb' ■ Let 

c = min^gg^) 0(a:;p,p'). Let T satisfy (1.5). 

Assume that D is admissible and dD is . If A.qb{p-iP') is finite and for all q e 

det{S,{E,{p,p')) - Sg{dD)) > 0, (1.14) 

then we have 

lim r^e^'"'''-esD,.eaB,«'eaB"^(^».2'') / (fv„-Wfg)dx 

T— foo J-R3\D 

/ N / ,N (1-15) 

TT / diamB \ j diamB \ 1 

" 2 ,eA,t(py) V2k-Ply ' Uk-P'lj ■ ^det{S,{E,{p,p'))-S,{dD)y 

Some remarks are in order. 

• The right-hand side of (1.15) is symmetric with respect to the replacement p ^ p' 
and p' p. This is a kind of reciprocity. 

• The quantity det {Sq{Ec{p,p')) — Sq{dD)) expresses some kind of information about 
the difference or deviation of the geometry between dD and Ec{p,p') at q & h.Qr,{p,p'). 

The following proposition says that Theorem 1.3 can cover convex obstacles. 

Proposition 1.1. (i) If D is convex, then D is admissible and Aqu{p,p') consists of a 
single point, (ii) Let q G Aqd{p,p') and assume that dD is contained in the half space 
{x — q) ■ Uq < 0, then (1-14) at q is satisfied. 

For the proof see Appendix. Thus as a corollary we obtain the following result. 

Corollary 1.1. Let B and B' satisfy [B U W] f] D = $. Let f = Xb and g = xb' ■ Let 

c = mmx(zdD(p(yX;p,p'). Let T satisfy (1.5). If D is convex and dD is , then (1-15) 
whose right-hand side consists of a single term is valid. 
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Note that, for the back-scattering case B — B', using the completely same argument 
as done in [12] in a bounded domain, we obtain 



hm rV^dist(9D,B) f ^^Adx 

T i-OO Jb 

71 ( diam B \ 



(1-16) 



E 



2\2ddD{p)J xedD,\x-p\=dgD{p) ^/Pdoil/d^Dip}^ 

where daDip) = inf^^eSD \x - p\; PdD{><; q) = {X - ki{q)) (A - k2{q)) and ki{q) and k2{q) 
denote the principle curvatures of dD at q with respect to Ug (see Appendix in [6]). Note 
that the Gauss curvature KQj:,{q) and mean curvature HgD{q) at q with respect to i/g are 
given by ki{q)k2{q) and {ki{q) + /c2(g))/2, respectively. 
The assumptions therein are 

• dD is C3; 

• T > 2dist {D,B); 

• the set of all points x e 5L> with |x — p| = dgoip) is finite and each point ^ in the 

set satisfies PQj:,{l/dQ£,{p);q) > 0. 

It is not assumed that D is admissible in (1.16) unlike (1.15). 

The quantity Paoi^/ dQD{p)]q) at g e dD with \q — p\ = ddoip) denotes a 'deflection' 
of the surface dD at q from the sphere |a; — p| = dgoip) since we know from, e.g.. 

Proposition 4.2 in this paper that Pdoiq', ^/doip)) = det {Sg{dBdjy(^p){p)) — Sg{dD)) , where 
Bdjj(p){p) = {x E 'R? \ \x — p\ < d]j{p)}. Thus formula (1.15) of Theorem 1.3 can be 
considered as an extension of (1.16) to the bistatic data case. See also Remark 5.2 for a 
comparison. 

After having Theorem 1.3, everyone wishes to extract the geometry of dD at all the 
first reflection points. The complete answer for general obstacle is not known, however, 
under the admissibility of D, one can obtain the following result. 

Theorem 1.4. Let q G Kqo{p^p') be known. Let T satisfy (1.5). Let B and B' satisfy 
\B\J B]r\ D = %. If D is admissible and dD is C^, then one can extract 

„ S,(dD)(A,(p) X A,{p')) ■ (A,(p) X A,{p')) 

and 

Kaoiq) 

from Uf on B'x]0, T[ with f = xb, where 

\q — x\ 



Note that Ag{p) x Ag{p') belongs to the tangent space of dD at q. For this see Lemma 
4.3 in Section 4. 

This theorem may suggest the following. 

• If one wishes to know the mean curvature at a first refiection point precisely, one 
should make the transmitter and the receiver approach as much as possible. It is because 
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Aq{p) X Aq{p') will disappear approximately at this time and thus the correction term in 
(1.17) can be ignored. 

• On the other hand, the Gauss curvature at the first reflection point can be extracted 
regardless of the position of a transmitter and a receiver at any time except for the 
condition [B UB'] nD = 0. 

As a corollary of Theorems 1.1, 1.2 and 1.4 we have the following result. 

Corollary 1.2. Assume that D is an open ball. Let T satisfy (1.8). Let B and B' satisfy 
U S ] n D = 0. Then, one can extract D itself from, Uf on B'x ]0, T[ with f — xb- 

The steps to reconstruct an unknown open ball D are as follows. 
Step 1. Determine c = minj;^^^ via Theorem 1.1. 

Step 2. Determine the unique point q in AgE,{p,p') together with i/q via Theorem 1.2. 
Step 3. Determine KQ£,{q) via Theorem 1.4. 



Then the radius and center of D are given by 1/ ^ KQu{q) and q — {1/ y KQD{q))i'q, 
respectively. 

The reconstruction problem of a spherical obstacle also has been considered in the 
frequency domain. For example, see [1] which employs a spherical wave as an incident 
wave and uses a low frequency limit for the reconstruction. 

The four steps described above give a constructive proof of a uniqueness theorem in 
an inverse obstacle problem in the sense that it does not make use of the uniqueness of 
the continuation of the solution of the governing equation of the wave. The following 
uniqueness result employs the bistatic data over a finite time interval and itself seems to 
be new. 

Theorem 1.5. Let Di and D2 be open balls. Let u-'^ be the solution of (1.1) with f — Xb 
and D = Dj. Let T satisfy (1.8). Let B and B' satisfy \B U b'\ nDj^^ for j = 1, 2. // 
u) = u} on B'x ]0, Tl then Di = D2. 

We refer the readers to [16, 17, 23, 24] for various uniqueness theorems for inverse 
obstacle problems for hyperbolic equations over a finite time interval. 

Another corollary from Theorem 1.4 is concerned with the determination of the direc- 
tions of principle curvatures at a point on dD. 

Assume that D is convex and dD is C^. Prom Proposition 1.1 we know that Agoip^p') 
consists of a single point. We denote the point by q{p,p'). We denote by p{9) and p'{9) the 
points rotated around the line directed Uq at q = q{p, p') counterclockwise with rotation 
angle 9 e [0, 27r[ of p and p'. Thus p(0) = p an p'(0) = p'. 

Then for all e e [0, 27i[wekaow thai A9D{pi9),p'{9)), Aq{p{9))-Aq{p'{9)), \Aq{p{9))x 
Aq{p'{9))\, (t>{q{9)]p{9),p'{9)) and VqOi q — q{p{9),p'{9)) are invariant with respect to 9. 

Let B{9) denote the open ball centered at p{9) with radius rj and B'{9) the open ball 
centered at p\9) with radius 77. We have \B{9) U s'(6')] r\D = provided [Bu'b'] CW = 
and D is convex. 

Then from (1.17) in Theorem 1.4 applied to / = f{9) and B = B{9) and B' = B'{9) 
we obtain the function of 9: 





where V(^) denotes the unit vector directed to Aq{p{9)) x Aq{p'{9)). 




Sq{dD){V{9))-V{9) 
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Now assume that Aq{p) x Aq{p') ^ 0. Then, V{9) attains all the tangent vector at 
q of dD and thus from the behaviour of HQ]:,{q;p{9),p'{9)) as a function of 9 one can 
determine all the directions of principle curvatures say, V^(^i) and ^(^2) with some 9i 
and ^2- Then we have 



H9D{q;p{9i),p'{9,)) + Hon{q;pi92),p'i9,)) _ \ 1 1 - Aq{p) ■ Aq{p') 




l + Aq{p).Aq{p')) 



HdD{q). 



Thus we obtain Hgciq)- 

Summing up, we have obtained the following result. 

Corollary 1.3. Let B and B' satisfy [5U 5 ] nU = 0. Let T satisfy (1.5). A ssume that 
D is convex and dD is ; q = q{p,p') is known; Aq{p) x Aq{p') ^ 0. Then, one can 
extract all the directions of principle curvatures, mean and Gauss curvatures, in other 
words, the shape operator at q of dD from Uf(^g) over B'{9)x ]0, T[ for all 9 e [0, 27r[, 
where f{9) denotes the characteristic function of B{9). 

A brief outline of this paper is as follows. Theorems 1.1 is proved in Sections 2 and 3. 
As described above, the key point of the proof is to derive (1.9) and (1-10) and those are 
proved in Sections 2 and 3, respectively. 

Theorem 1.2 is proved in Subsection 5.1. The proof contains an explicit characteriza- 
tion of the first reflector in terms of the bistatic data. See Remark 5.1 for the resulted 
procedure to determine all the first reflection points. 

Theorem 1.3 is proved in Section 4. The key point in the proof of (1.15) as well as 
(1.16) is to identify the term which contains the leading term of the indicator function. 
See (4.1) for the term. We found that the one of two reflection arguments developed in [18] 
works for the purpose. It is based on the reflection across dD and a pointwise estimate 
of e° near dD, that is the use of the admissibihty of D. The argument is presented in 
the proof of Lemma 4.2 in Section 4. Note that another reflection argument used in the 
proof of (1.16) is free from the admissibility assumption, however, can not be applied to 
the case when f ^ g. 

Theorem 1.4 is proved in Subsection 5.2. The proof is based on an asymptotic formula 
which is a consequence of Theorem 1.3 and an explicit formula of the determinant of 
the difference of two shape operators at g e Aqd{p,p') as derived in Subsection 7.3 of 
Appendix. 

In the final section we give a conclusion of this paper and comments on further prob- 
lems. 



2 An upper bound of the indicator function 

Define 

J{r; f:9)^ [ i^Vf ■ Vvg + T%fVg)dx. (2.1) 
We have the following expression of the indicator function. 
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Proposition 2.1. It follows that 



,-Tr 



Proof. Prom (1.2) and (1.4) we have 

/ ifva — Wfq)dx — I '^^^vAS -\- e~'^'^ [ Ff(x,T)vodx. 

Rewrite 

JdD OV JdD OV JdD OV JdD OV 

Integration by parts yields 

f dS=J{T-J,g). 
JdD du 

On tlie other hand, (1-13) yields 

Furthermore it follows from (1.2) and (1-11) that 

dZf 

Ir3\d ' 



P dZ f 



Now from this together with (2.3)-(2.6) we obtain (2.2). 
□ 



Lemma 2.1. As r — > oo 



Proof. From (1.2), we obtain 
Since 



Lj^^^r^v, + T^ZfV,)dx 



2t2 



1l, ,2 \9\' 



4^4' 



from this we obtain 
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and thus 



Similarly it follows from (1.11) that 

/ _{\^Zf\^ + T^\Zf\^)dx< \ f 



(2.8) 



\^dx. 



(2.9) 



A combination of (2.8) and (2.9) and the estimate ||-^/||l2(j^3\;d) = 0{t) yields (2.7). 

□ 



Thus a combination of (2.2) and (2.7) gives 

/ , _(M - njjg)dx = J{t- f.g)+ I _(VeJ ■ VeJ + Thyg)dx + 0\ 



For the second term in this right-hand side we have the following estimate. 
Lemma 2.2. As r 



T-^e-^^). (2.10) 



oo 



(2.11) 



Proof. It is an application of the trace theorem twice and integration by parts. More 
precisely, choose v e H^{R^ \ D) in such a way that 



/2 



(2.12) 



where C > and is independent oivf. Integration by parts (or the weak formulation of 
(1.13)) yields 



/ -^VfdS=- 

JdD OV 

A combination (2.12) and (2.13) gives 



(Ve*] • Vv + T^e^fv)dx. 



(2.13) 



-K^VfdS 

dD OV 



< C(||Ve}||^2(j^3\:D) +'^14llL2(R3\B))lk/l9ollHV2(aD)- (2.14) 



Since 



max 



, -^VfdS 

dD OU 



|V6?|L.(^3\^),r||4||,.(^3\^)) < (/^^ _(|V6jr + r^|4ncix) 
it follows from (2.14) that 

<C{l + r) (/^3_(|V6jp + r2|6jr)dxy^'||^^|,^|U: 
From this together with (2.6) for f = g, we obtain 



1/2 



(dD)- 



\dDy 



(2.15) 
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By the trace theorem, we have 

<c^'(liv^/lli.(z.) + ll^/lli.(z.)), 

where C" > is independent of Vf. Now from this together with (2.15) and the trivial 
estimates 

max(||V^;/||i2p),T2||^;/||i2p)) < J{T;fJ), 

yields (2.11). 

□ 

Therefore (1.9) with jii = 2 follows from (2.10) and (2.11) together with the following 

estimate. 

Lemma 2.3. We have, asT — > oo, 

J{T;f,g) = C>(re-^"^^^xe9o,.e8B,.'eaB'^(=^;2/>s/')). (2.I6) 

Proof. It follows from (1.3) and (2.5) that 

J(r-J,g)=(^)' [ dsj K{x,y,y')dydy', (2.17) 

\47r/ JdD JBxB' 



where 

kr{x,y,y') = + r] JK^^hj^e'^(\^-y\M^-y'\), (a:,y,y') edDxBxB'. 

\\x — y\ J \x — yl'^lx — y'\ 

(2.18) 



Since we have 



inf d>(x: y, y') — min (h(x: y, y') 

xedD,yeB,y'eB'^^ '^'^^ xedD,yedB,y'edB'^^ '^'^^ 

and BnD = B'nD = 0, from (2.17) and (2.18) we obtain (2.16). 
□ 



3 A lower bound of the indicator function 

3.1 A reduction to a convex obstacle and the proof of (1.10). 

Rewriting the second term in the right-hand side of (2.10) with (2.6), one has 

, _Jv„ - Wfq)dx = J[t\ j,q) + I 



/ (fvg-Wfg)dx = J{T;f,g)+ / -^v,dS + 0{T-'e-^^). (3.1) 

JR^\D JdD OV 



In general, we do not know the signature of the second term of the right-hand side of (3.1), 
however, we know that the function under integral is nonnegative at a special point on 
dD by virtue of the following lemma which is an application of the maximum principle for 
differential operator A — and a reflection argument in [18]. It corresponds to Lemma 
3.7 in [18] in which Vf in (1.13) is replaced with —e~'^^''^ for a a; e -S"^. 
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Lemma 3.1. Let q e dD be a point of support of D, i.e., D is contained in the half-space 
X ■ Ug < q ■ Uq. We have 

^(5) > 0- (3.2) 

Proof. First, we prove that, for all x e \ D, 

e){x) > -Vfix). (3.3) 
Define wj ^ + Vf E H\R^ \D). It holds that 

{A-r^)wj = -f mR^\D, 
= on dD. 

Since (A — T'^)w'j < in \ D, from the weak maximum principle for operators of 
divergence form (Theorem 8.1 in [6]), we obtain 

inf w°f > min (inf w% 0). (3.4) 

where Bji denotes an arbitrary open ball centered at the origin such that DUsupp / C Bji 
and Sr = dBji. Since 6° decays as — > oo uniformly with respect to wc have 

inf^^ wj — > as R — > oo, and thus, from (3.4) we obtain w^{x) > for all x Ell^\D. 
This completes the proof of (3.3). 

The equality in (3.3) holds for x — q. This implies the following inequality for the 
normal derivatives: 

f (3.) 
Second, we prove that, for all points that satisfy x • Uq > q ■ Uq, 

6?(x) > -Vf{x'), (3.6) 

where x' is the image of x under reflection across the plane x • Uq — q ■ Uq. Since x — x' 
on the plane x ■ Ug = q ■ rig, (3.3) shows that (3.6) is satisfied there. Note that also 
v'j^{x) = Vf{x') satisfies (A — t'^)v'j = —f{x') < 0. Applying the weak maximum principle 
to {wjy = e^f + v'f in the half-space x ■ rig > q ■ rig, one obtains as before (3.6) holds 
throughout the half-space. This completes the proof of (3.6). 

Since the equahty in (3.6) holds for x — q,\t follows as before that 

fM>-|;WM}U,-|fte). (3.7) 

Now a combination of (3.5) and (3.7) yields (3.2). 
□ 

The following lemma is an easy consequence of the C^-regularity of dD and thus the 
proof is omitted. 

Lemma 3.2. Let q e Aqd(p,p'). Then, there exists an open ball D contained in D such 
that q e Ag^(p,p') and thus min^^^^ 0(x;p,p') = mm^^gD(l)(x;p,p'). 
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Let u — Uf denote the weak solution of the following initial boundary value problem: 

d^u - A-S = in (R=^ \ ^) x ]0, T[, 
u{x,0) = OmR^\D, 
dtuix,0) = f{x) mR^\D, 
u^O on9^x]0, T[. 



Define 



Wf{x, t) = J e'^'^u{x, t)dt, X eR^\D, T > 0. 



Lemma 3.3. We have 

/^3^_(M - y^f9)dx > l^,^-gifv9 - ^f9)dx + 0{T-'e-^^). (3.9) 

Proof. Let Zf e H\R^\D) solve 

(A - T^)Zf = Ff{x, t) in \ 
= on9L>, 

where 

Ff{x, t) = dtUfix, T) + TUf{x, T), x G \ Z). 

Similar to Zf which is the solution of (2.2), we have 11-^/11^2(^3^5) ~ And similar 

to (1.12) for Wf, Wf has the form 



where satisfies 



Wf^Vf + ^f + e-^^Zf, 

{A-T^)^} = QmR^\D, 
= —Vf ondD. 

Thus, we have 

Wf-Wf^ (e? - ej) + e-^^{Zf - Zf) in R^ \ D. 

Since Vf>Ooia dD and ^f{x) — > as \x\ — > oo, by the maximum principle for the 

modified Hclmholtz equation in R^ \ D, we have — < Vf in R'^ \ D and thus — < — ej 
on dD. Again by the maximum principle for the modified Helmholtz equation in R^ \ D, 
we obtain — e° < — e° in R^ \ D. Therefore we obtain 

Wf-Wf>e-^^{Zf-Zf)mR^\D. (3.10) 
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Since both supp g and supp / are contained in R^\Z^ and thus combining this with (3.10), 
we obtain 

L.y,^^^3 - Wfg)dx - Ijyfvg - Wfg)dx > e"^^ L^-i^f - Zf)gdx. 
Jb?\d jr3\d Jr3\d 

From the L^-bounds for Zf and Zf, we see that this right-hand side has the bound 

0(r-ie--^). 
□ 

Since D is convex, every point q e dD is a point of support of D and thus, from (3.1), 
(3.2) and (3.9) we obtain 

/^3^_(M - Wf9)dx > J{T;f,g) + 0(r-V--^), (3.11) 

where 



J {r; /, g) = /- (Vv/ • Vvg + T^VfVg)dx. 
Jd 



Now everything is reduced to give a lower estimate for J{T;f,g) as r — > oo. For this 
and the future use of it in the sound-hard obstacle case we give the estimate for J(t; /, g) 

for general D. 

In the following lemma we do not assume that D is convex. 
Lemma 3.4. There exist positive constants C , fi and tq such that, for all r >Tq, 

min^gac -Pix;P,p') ^-rivW) J (^^. > C. (3.12) 

We give the proof of this lemma in the next subsection. 

It follows from (3.11) and (3.12) lor D — D that there exist positive constants C and 
Tq > such that, for all r > Tq, 

^2+Mgrmin,gao,^gaa,^,gaa,<^(r.;2/y) J^^_{fVg - Wfg)dx > C 

provided T satisfies (1.5). This completes the proof of (1.10). 
3.2 Proof of Lemma 3.4. 

In this subsection we never assume that D is convex. Let A{x, t) be an arbitrary positive 
function of a; G dD with parameter r > 0. For another function i?(x, r), in the following, 
B{x, r) = 0{A{x, r)) as r — )■ oo and uniformly with respect to a; G dD means that there 
exist positive constants tq and C independent of x G dD such that, for all r > tq and 
X G dD we have \B{x,t)\ < CA{x,t). 

The proof of Lemma 3.4 starts with having the following expression. 

Lemma 3.5. There exists a positive constant C such that, as r — > oo 

J(T\f,9) 

iT-^ \ T J \ T J JdD \X — prlX — p'\ \ T\X — P\ 



At^ V T J \ T J JdD \x — pp|a; — p'\ \ t\x — p\ 

_|_Q^^-lg-T infa,gao{ds {x)+dgi {x)) (^^-Ct rain^^gu ds (x) _|_ g-Cr min^^go dg, (x) ^ 



(3.13) 
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Proof. By [12], we have, as r — > oo and uniformly with respect to x e dD, 



2r 

and 



where C is a positive constant. 



p-TdB{x) / 1\ / 1 

= — — -(77-- 1 + 



p-TdBix)^yl+{^r,/dB{x)) / / 2^ 



and 



( J+),(a;, r) = V - - + 7\ ■ 

t\x — p'\ V T/ T^|X— p'l 

Prom these and (2.5) we obtain 

Jir;f,9) = ^ / (J+)/(a:,r)^^^^(J+),(a:,r)d5. 

j^Q j'^-lg-r inf^gao (ds {x)+dgf {x)) (^^-Ct min^gao da {x) _|_ g-Cr min^ggo (^s' i^) ^ j 

This yields (3.13) since we have as r — )■ 00 and uniformly with respect to x e dD, 

p-T[dB{x)+dg,{x)) / IX/ 1 \ / 1 \ 

WU^,r)iJi),i.,r) ^ ^,|^_^||^_^| (, - -) (V - -) (1 + ^) 

j^Q{j-^Q-r{dB{x)+dBi{x))^^-TdBi(x)C _^ g-Tds(a;)C-j-j _^ Q^^-4g-T(ds(a;)+ds, (a;))(l+C) -j _ 



jaD \x — v\ \x — v \ 



□ 

Now we give a lower estimate of J(t; f,g) as r — V 00 by using (3.13). 
Define 

IdD \X — pY^Ix — p 

where m = 2, 3. 

Since dsix) = \x — p\ — rj for x & \ B and dB'{x) — \x —p'\ — r]' for a; e \ B' , it 
follows from (3.13) that 

e-<-V(.; /, ,) ^ ^ (, - 1) (,' - 1) (/.(.) + l/3(x)) ^^^^^ 

_|_Q^^-lg--rmina,gao</>(a:;j'y)^g-CTmin^gao<is(a;) _j_ g-Cr min^gao dg/ (a;)) 

Since 

^ [ {p-x) \ (m-3) ^ (p - x) • (p^ - x) 

] |x — — p'l J \x — p\^\x — p'\ \x — p\^\x — p'\^ 
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and 



we have 



{p — x) ■ V(p{x;p,p') 

^-T<j>(x;p,p') 



\x — p\ + 



ip — x) ■ {p' — x) 



\x — p' 



_3) / , ^'^ '7"' „rfx+ / (P-^)- (p--^) -.,(.„.,^ 

Jd \x — py^\x — pf] Jd \x — py^\x — p'\'^ 



7. 



(m-3) 



X — p\"''\x — p 

T(j){x;p,p') 

-dx + 



D \x — p\'^\x — p'\ J D \x — pY^\x — p' 



X 



+^ / I ] — ^ Trdx + r / f— 



ip- x)- ip' -x) ^ 
p\"^\x — p'p 



r4>(x;p,p') 



This yields 



/2(t) + ^/3(t)=t|^|i + 



(p-x)- {p' - x) \ e-^^(^;P'P') 



+ 



X — p\\x — p'\ J |a;— p||a;— p' 

1 



I 1 + 

D \\X — p\ \X—p'\l \X — p\^\X — p 



/1 2 



(3.15) 



+ 



1 I (p^^^Hy^g-r^(.;py)^^_ 
T JD \x — p\^\x — p'P 



Lemma 3.6. Let D he an arbitrary nonempty hounded open set. If p and p' he arbitrary 
points in'R?\D such that n D = 0, then 



Cn{p,p')=iTii{l + 



{p — x) ■ [p' — x) 



x^D y \p — x\\pf — x\ 



> 0. 



(3.16) 



Proof. It is easy to see that < Cd{p,p') < 2. Assume that Cd{p,p') — 0. Since we have 
the identity 

2 



1 + 



(p — x) ■ ip' — x) 1 



p — X ^ p — X 



\p — x\ \p' — x\ 



, xj^p,p', 



there exists a sequence {xn} in D such that 



P-Xn , P - Xn 



\p-Xn\ \p'-Xr, 



0. 



(3.17) 



Since D is compact, choosing a subsequence of if necessary, one may assume that 
{xn} converges to a point y G D. Since p,p' ^ y by assumption, it foUows from (3.17) 
that 

p-y , p' -y 



+ 



p-y\ b' - y\ 



0. 
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This gives y G \p,p'] and thus \p,p'] H D ^ 0. This is a contradiction. 
□ 

A combination of (3.15) and (3.16) gives 



Hr) + -h{T) > tK{t) I e-^'^^'^'P'P'Ux, (3.18) 

T JD 



where 



ddD{p)ddD{p') \daD{p) daoip') j ddD{p)ddD{p') r dduipfdQDip'Y r"^' 
Lemma 3.7. Letp,p' e R^. Then, there exists a number /i such that 

hminfr^e^^^'^-eaD^Co^ipy) /■ e'^'^^'^'P'P'Ux > 0. (3.19) 
T — >-oo Jd 

Proof. Let xq G dD be a point such that (f){xo;p,p') = miiixedD (l){x;p,p')- Since |p — x| < 
\p — xo\ + \xq — x\ and \p' — x\ < \p' — xo\ + \xo — x\, we have 

(f){x;p,p') < (f){xo;p,p') + 2\xo — x\, Vx G R^. 

This gives 

Jd ~ Jd 

In [15] we have aheady known that 

hminfr^ / e'^^^^'-'^'Ux > 0. 



T >00 



D 



Thus (3.19) is vahd for = 3. 
□ 

Now it follows from (3.14), (3.18) and (3.19) that (3.12) is valid. 

4 Asymptotic behaviour of the indicator function 

First we claim that 

/ 3,-(M - ^f9)dx = 2 J(r; /, ^)(1 + ©(r'V^)) + ©(r-^e"^^). (4.1) 

This is a consequence of the following asymptotic formula and (2.10). 
Proposition 4.1. As r — > oo, 

/^3^_(V6j • Vel + T^el)dx = J{t- f, g){l + 0{t-'/')). 
Proposition 4.1 is a direct consequence of the following lemmas. 
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Lemma 4.1. Let D be an arbitrary nonempty bounded open set. If B and B' be arbitrary 
open balls such that \B \J B] r[ D — 0, then there exist positive constants C and tq 
independent of r such that, for all r > Tq, J{t; f,g)>0 and 

Jir]f,g) > Ct^ I dx I e-^'^^'^'y^y'^dydy'. 

JD JbxB' 

Proof. From (1.3) and (2.1) we have 

V47r/ Jd JbxB' \x — y\\x — y'\ 



-T<p{x;y,y') 



(4.2) 



where 



Krix,y,y') = 



1 ^ \ {y-x) _ {y' - x) ^ ^2 



\\(x — y\ J \\x — y'\ J \x — y\ \x — y'\ 
Since BOD — B'dD — $, there exist positive constants Ci and C2 such that 

K^ix, y, y') > f 1 + • - C^T - C2, V(x, y,y') e D x B x B'. (4.3) 

V \x-y\ \x-y'\J 

Here we claim the following estimate: 



Co{B, B') ^ ^ inf f 1 + • ^5-4) > 0. 

{x,y,y')eDxBxB' \ \y — x\ \y' — x\ J 



(4.4) 



This is proved as follows. 

It is easy to see that < Cb{B, B') < 2. Assume that Cd{B, B') — 0. Since we have 
the identity 



1 + 



y — X y' — X 1 

ll/ — ,t| \y' — x\ 2 



y — X ^ y' — X 



\y — x\ \y' — x\ 



, X ^ y,y', 



there exist sequences {xn} in D, {yn} in B and {y!^} in B' such that 



Vn Xji ^ y^ Xji 



I Vn Xji I I Xji I 



0. 



(4.5) 



Since D x B x b' is compact, choosing a subsequence of {(a;„, |/„, |/^)} if necessary, one 
may assume that {{xn,yn,y'n)} converges to a point (a;*, |/*, 7/!|,) G D x B x B . Since 
y*iy'<^ 7^ X* by assumption, it follows from (4.5) that 



y* Xjf y^ x^ 



\y* xj\ ||/^ X, 



0. 



This gives e {sy^ + (1 - «)?/* | < s < 1} and thus [S U S ] n D 7^ 0. This is a 
contradiction. 

Thus, applying (4.4) to the right-hand side of (4.3), we obtain, for a sufficiently large 
To > 



r 



inf i^.(a;,y,|/') > -Cn{B,B'), Vr > r^. 

{x,y,y')eDxBxB' Z 
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A combination of this and (4.2) ensures the vahdity of Lemma 4.1 with 

t^Cd{B,B') 



c 



4:7tJ 2dist {D, B) dist (D, B') ' 



□ 



Lemma 4.2. Assume that D is admissible and its boundary is C^. Then, there exist 
positive constants C and Tq independent of t such that, for all r > Tq, 

I -(Ve". • Ve° + T^yMx - J(t; f, g) < Cr'^l'^ [ dx [ e-^^^'^'^'^'^ydy'. 
Jn3\D J y J y JbxB' 

Proof. Let < 6 < 60- Let = 05 be a smooth cut-off function, < 0(x) < 1, and 
such that: (p{x) = 1 if dooix) < S and (p{x) = if daoix) > 26; \V(j){x)\ < CS~^; 
|V20(a;)| < C5-\ 

For * — f,g define (e°)^(a;) = (f){x)e1{x'^) for x e D and vl{x) — (l){x)v^{x'^) for 
xe R^\D. 

The Lax-PhiUips reflection argument starts with the following expression: 

/r3\b^^^/ ■ + ^'44)dx = At; f, g) + J^^_ eJlA - r'^dx. (4.6) 

In the proof the following relationship between vl and v^^ and the boundary condition for 
e° are essential: 

dvl dv^ 

vl — V* — — e° ondD. 

Another device is the following differential identity which is a consequence of (4.15) 
in [18] (see also Appendix 1 in [12]): 

(A - r'){vl) = <j^{x)Y.a^Axmd,v,){x^) + E + (A0)(x)t;,(xO. 

i,j 3 

where ajj(x), i,j = 1, 2, 3 are in a neighbourhood of SD, independent of and Vg and 
satisfy 

3C> OVx e (t){x)\aij{x)\ < Cdaoix)] (4.7) 

each hj{x) has the form 

H^) = E hk{x)dk(l>{x) + dj{x)(t){x) (4.8) 

jk 

with bjk{x) and dj{x) which are and C° in a neighbourhood of dD, respectively and 
independent of and Vg. 
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This together with the change of variables x = y'^ yields 

j e]{x){A-r^)vl{x)dx 



D 



y i,3 j ) 

(4.9) 

where J{y) denotes the Jacobian. 

Hereafter, we give an estimation for each term in the right-hand side of (4.9) point- 
wisely, without making use of integration by parts further. This idea is exactly same as 
the proof of Lemma 3.3 in [18]. This is different from the back-scattering case, see also 
the proof of Lemma 4.2 in [18] and Appendix 1 in [12] for the comparison. 

Since D is admissible we have, for all y G D with ddEi{y) < S' and all r > Tq 

|e5(/)| <C [ e-^l^-^ldx, (4.10) 

where C is independent of y and r. 
It is easy to see that, for all y E D, 



T-^\didjVg{y)\+T-'\d,v,{y)\ + \vg{y)\ < C j ^e-^\y-^'kx' . 
Prom this, (4.10), (4.7) and (4.8) and the choice of (f) we obtain, for all y E D, 

^'M) \W)i:<^iAfMd,v,){y) + Y.Wmv,){y) + {A^){f)v,{y) 

y i,3 3 

JbxB' 



(4.11) 



Choosing S — t""^/^, we have Sr^ + + = 0{t^^^) as r — > oo. Now from (4.9) 
and (4.11) we obtain the desired conclusion of Lemma 4.2. 

□ 

Thus everything is reduced to studying the asymptotic behaviour of J(r; /, g) as r — > 
oo. For this purpose we employ the asymptotic formula (3.13). Note that in Section 3 we 
made use of the formula to give a lower estimate of J(r; /, g). Here using the formula, we 
determine its leading term as r — > oo. 

From (3.13) we see that the asymptotic behaviour of the following integral is the key: 

/ , • , ( 1 + e-^(--')d5.. 

JdD \X — p\^\x — p\ \ T\x—p\J 

Proposition 4.2. Assume that Aqd{p,p') is finite and 

det{Sg{E,ip,p')) - S.idD)) >Oyqe A9Dip,p'). 
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Then, we have 

li^ ^^r4>(,;p,p') f - ( 1 + —1—] e-^t>^^'^'^'^S, 

— ^oo JdD \X — p\^\x — p'\ \ T\x—p\J 

TT 1 



= V - 

,eAeo(j>,P') I? -P'l ^Jdet{S,{E,{p,p')) - Sg{dD)) 

Theorem 1.3 is a direct consequence of this together with (3.13) and (4.1). 
In the following subsection, we describe the proof of Proposition 4.2. 

4.1 Proof of Proposition 4.2. 

We employ the Laplace method and so one has to compute the Hessian of the real phase 
function dD 3 x i — > (j){x;p,p') at all the points on dD where it takes the minimum value. 

Let q e Aqd{p,p'). One can choose a local coordinates system a = {ai,a2) around q 
on dD in such a way that x e dD around q has the form 

X = Xq(a) ^ q + (7iei + 0-262 + f{cr)vq, \a\ < 5, 

where 5 is a sufficiently small positive number independent of 5; Ci and 62 are two unit 

tangent vectors at q to dD which are perpendicular to each other and ei x €2 — Vq'-, 
f ^ fqe C3(R3) and satisfies /(O) = 0, V/(0) = 0; takes the form 

-fiei - 7262 + 



where /i = df/dai and /2 = df/da2; dS = ^jlVW^rJida . 
We have, for z — p,p' and j = 1, 2, 

d 1 

'Xq{a) -Z\ = —J- {{Xq{a) - Z) ■ Bj + fj{a){Xq{a) - Z) ■ Uq] (4.12) 



dUj \^q{'^) ~ ^ 

and thus 

d^ 



\xq{a) - z\ 



dakdaj 

^ {{Xqjcr) -z)-ek + fk{cr){Xq{a) - z) ■ Uq}{{Xq{a) - Z) ■ Bj + fj{a){xq{a) - z) ■ Uq} 

\xq{a) - z\^ 

This yields 

^^^^0(xq((7);p,p')U=o = K(l'^P,p')Skj - akj(q;p,p'), (4.13) 
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where 

Kt,P,p) = 1 r + T-, r, 

\p - (l\ \p - q\ 

1 1 (9V 
akj{q;p,p') = 1 ;A ■ e^A ■ Cj + — -A! ■ e^A' ■ - (0) {A-v^ + A!- Uq) 

\p-q\ \p-q\ dukda-j 

and 

Note that \{q;p,p') = \{q;p',p) and akj{q;p,p') = akj{q;p',p). 

The foUowing lemma corresponds to SneU's law in geometrical optics. 

Lemma 4.3. Let p and p' be arbitrary points in R'^ \ dD such that [p,p'] fl dD = 
Let q G Aqd{p,p'). Then, we have: (i) the vector A + A' is parallel to Uq; (ii) A ■ Uq = 
A! ■ 0. In addition, if p and p' are in \ D, then we have 



V, 



1 



^2(1 + A-A') 



[A + A'), (4.14) 



that is, the unit inward normal to Ec(p,p') at q with c — (l){q;p,p') (see (A. 17) in Ap- 
pendix) coincides with the unit outward normal to dD at the same point. 

Proof. Since the function a i-> (f){xq{a);p,p') takes its minimum at o" = 0, it follows from 
(4.12) that 

A • ei + A' • ei = 

(4.15) 

A • 62 + A' • 62 = 0. 

Write 

A = aei + /3e2 + ^Vq, 

(4.16) 

A' = a'ei + j3'e2 + -I'uq. 

(4.15) is equivalent to a + a' = and /9 + /?' = 0. Then we have 7^ = 7'^, that is, 7 = 7' 
or 7 = —7'. Assume that 7 = —7'. Then from (4.16) we have A + A' = 0, that is, 

q-p ^ q-p' ^ Q 



\q — p\ \Q — P 

This means that q G [p,p'] and contradicts the condition [p,p'] fl dD — 0. Thus 7 = 7' 
and we have 

A + A' = 271/, (4.17) 

and 

^Vl + A- A' 

Since | A + A'p = 2(1 + A ■ A'), from the argument above we have 1 + A • A' > 0. Thus 
from (4.17), one gets A ■ Uq ^ and A' ■ Uq ^ 0. 

Note that if p and p' are in R^ \ D, 7 has to be negative. The reason is the following. 
Assume that 7 > 0. Then — (A + A') is directed to —Uq. Since dD is C^, one can find a 
sufficiently small s > such that q' = q — s{A + A') e D. Since p G R^\D, one can find a 
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point p" e dD on the segment with endpoints q' and p. Note that —{A + A') is directed to 
the unit inward normal to Ec{p,p'). This together with the condition c > \p — p'\ ensures 
that both q' and p are in the domain enclosed by Ec{p,p') and thus by the convexity of 
the domain, we have 4>{p"',p,p') < c. However, since p" G dD, we have (f){p";p,p') > c. 
This is a contradiction. Therefore one gets 7 < and now it is easy to see that all the 
conclusions are valid. 
□ 

It follows from (i) in Lemma 4.3 that 

* G jEj * — * ^ fc * 



Thus, we have 



,^4. * G^-^i * 



\p-q 

^ Tj- — r I ^ • ^kA ■ Sj - — -A ■ BkA ■ Cj 



\p — q\ \p' — q\J \p' — Ql 

= X{q;p,p')A ■ BkA ■ Cj - ^ ^^ A ■ e^A' ■ Cj 



and hence 



- — - — :A ■ BkA ■ Cj + — rA' ■ BkA' ■ Bj = X{q;p,p')A ■ BkA ■ Bj. 
\p — q\ |P ?| 

Changing the role of p and p', we also have 

I ^ I A • BkA ■ Bj + — rA' ■ BkA' ■ Bj = X{q; p,p')A' ■ BkA' ■ Bj. 
\p — q\ \p — Q\ 

Prom these and (4.14), we obtain another expression for ajk{q '■ p,p')' 

ajkiq;p,p') - (A • BkA ■ e, + A' ■ BkA' ■ e,) + ^2{1 + A ■ A')-^{0). 

This together with (4.13) implies that 

Q2 



X 



^^^^^.0(^,(<t);p,p')U=o = ^|'2{1 + A■A') 

, ^^^'^'^'^ - \{A ■ BkA . B, + A' • BkA . Bj)\ - -P-{^) 

^2(1 + A • AV ^ ^ 9(7fe(7/ ^ 



(4.18) 



From (A. 2) we have 



T 
\ ^2 



S,{E,{p,p'))(b^ 62) 



, ^^^'^'^^^ , (S,k -liA-BkA. Bj + A' . BkA . Bj)) 

^2(1 + A • A') ^ 2^ "Jj 



j^l,2;fc-^l,2 
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and we know 



(el 



V el 



Sg{dD) ( ei 62 ) 



da j dak 



(0) 



(4.19) 



Thus we obtain the following formula which gives the geometrical meaning of the Hessian 

of (j){xq{a);p,p') at a = 0. 

Lemma 4.4. Let p and p' be inIl^\D such that \p,p'] fl dD = 0. Let q G Aqd{p,p') and 
c = 4'{(1]P)P')- Then, c> \p — p'\ and we have 

(el 



Vl<P{x,{ay,p,p%=o = ^J2{1 + A■A') {S,{E,ip,p')) - S^idD)) { e, 62 ) 



T 

V 62 



(4.20) 

Since (f){xq{a);p,p') takes the minimum at o" = 0, from (4.20) one concludes that, for 
all tangent vectors v ai q 



{Sg{E,{p,p'))-Sg{dD))vv>0. 



(4.21) 



Thus if (1.10) is satisfied, then from (4.21) one knows that Sq{Ec{p,p')) — Sq{dD) is 
positive definite on the tangent space at q and from (4.20) 

det V^0(x,(a);p,p')U=o = 2(1 + A • A')det {Sq{E,{p,p')) - Sq{dD)) > 0. 

And also from (4.14) we have 

k-pl 



{p-q)-Vq^ 



^/2 



l + A-A!. 



Now Proposition 4.2 is a direct consequence of the Laplace method [2]. 



5 Proof of Theorems 1.2 and 1.4 

5.1 Extracting the first reflector: proof of Theorem 1.2 

Let c > \p — p'\. Theorem 1.2 is based on the following proposition which gives a charac- 
terization of a first refiection point q between p and p' in terms of the minimum length of 
the broken path connecting pio q and q to p' -\- s{q — p') /\q — p'\ with a fixed small s > 0. 

Proposition 5.1. Fix < s < rj'. Let c — miUxedD p, p') . We have: 
(i) if p' + s{uj;p,p',c)uj belongs to dD, then 

min 4>(x;p,p' + scu) = c — s; 

xedD 

(a) if p' + s{uj;p,p', c)uj does not belong to dD, then 

min (f){x]p,p' + suj) > c — s. 
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Thus, one has the following characterization of the first reflector: 

^doiPiP') = {p' + s(uj;p,p', c)uj I min Mx^p^p' + su) = c — s, we S"^}. 

xedD 

Proof. Set p"{u!) = p' + scu. Let x e dD. We have 

(l)(x;p,p"(uj)) ^\p-x\ + \x-p"(uj)\ 

= \p-x\ + \{x-p') + {p' -p"{uj))\ (5.1) 
> |p — x| + |x — p'\ — \p' — p"{uj)\ — (f){x\p,p') — s. 

This gives 

min j9, > c — s. (5.2) 

x^dD 

Now we describe the proof of (i). Noting that s < s{u!;p,p',c) and p' + s{uj;p,p',c)u! e 
Ec{p,p'), we have 

(/)(p' + s{uj;p,p', c)uj;p,p"{uj)) 

= \p-{p' + s{uj;p,p',c)u) \ + \ {p' + s{uj;p,p',c)uj)-p'{uj)\ 

= |p - (p' + s(uj;p, p', c)uj) I + (s(uj;p, p', c) - s) 

= \p-{p' + s{uj;p,p',c)uj) \ + \ {p' + s{uj;p,p',c)uj) -p'\ -s 

= (/)(p' + s{uj;p,p', c)uj);p,p') -s^c-s. 
Thus if p' + s{uj;p,p', c)u! belongs to dD, then the inequahty in (5.2) 

min 0(x;p,p"(a;)) >c-s 

x€dD 

never occurs. Thus it must hold that mm^^QD 4>{x;p,p"{u;)) = c — s. 

Since we have always (5.2), the statement of (ii) is equivalent to the one that: if 
minj;^^^) 0(x;p,p"(a;)) —c — s, then p' + s{uj;p,p',c)uj e dD. 

Assume that miUx^dD 4'{x;p,p"{uj)) = c — s. Choose a point x £ dD such that 
(j){x;p,p"{uj)) = c—s. Then, from (5.1) we have c > (f){x;p,p'). Since c = mina^gau (f){x]p,p'), 
from this we obtain 4>{x;p,p') = c. Thus we have 4>{x;p,p"{uj)) = 4>{x;p,p') — s. This is 
equivalent to 

\x-p"{uj)\^\x-p'\-s. (5.3) 

Since x is outside the open ball B" centered at p' with radius s, one can find the unique 
point x' on dB" such that |a; — x'l = miny^gB" \y — x\. Since we have |x— p'| = — x'| +s, 
from (5.3) we obtain \x—p"{u) \ = \x—x'\. Since ^"(0;) G dB", it must hold that p" (a;) = x' 
and thus x — p' + s{uj;p,p', c)uj. This gives p' + s{uj;p,p', c)uj e dD. 
□ 

Now we are ready to describe the proof of Theorem 1.2. 

In what follows we denote the open ball centered at a point z and with radius p by 
Bp{z). Since B^i_s{p' + su) is contained in B' , from Uf on S'x]0, T[ with / = xb, one 
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gets Uf on Bn'-s{p' + su!)x ]0, T[. By Theorem 1.1 we obtain mfxeOD 4>{x;p,p' + sou) for 
each a; e S'^. Thus, from Proposition 5.1 we obtain AgD{p,p') itself. Prom formula (4.14) 
one gets Ug at given q G AgD{p,p'). 

Thus one can complectly determine the first rcfiectors between p and p' using the 
bistatic data Uf on B'x ]0, T[ for f = Xb and sufficiently large and fixed T. In particular, 
note that B is /ixec?. 

Remcirk 5.1. We summarize how to detect the points in h.QD{p,p') C Ec{p,p'). 
Step 1. Compute wj on B' with f — xb from the data Uf on S'x ]0, T[. 
Step 2. Pix s with < s < i]'. 
Step 3. Choose a direction u E S"^. 

Step 4. Compute the following integral with g = xb^,_^{p'+su>)'- 



/ Vgdx — / Wfdx. 



lB^,_^{p'+suj) 

Step 5. Compute the following quantity: 



— lim — log ( / Vndx — [ Wfdx] — n — n' + s. 

r-^oor \Jb JB^,_^ip'+su;) J 

Step 6. If the computed quantity in Step 5 is equal to c — s, then p' + s{iu;p,p', c)u G 
Aqjj{p,p'). If not so, then choose another cu and go to Step 4. 



5.2 Extracting the geometry of dD at a first reflection point: 
proof of Theorem 1.4 

In this subsection, we consider the case when a point q G Aqd{p,p') is known. We use the 
notation A and A' instead of Ag{p) and Aq{p'), respectively for simplicity of description. 
The aim of this subsection is to extract the geometry of dD at q. The main idea is to 
replace p' in AdD{p,p') with p' + sA' with a small s > 0. The advantage of this idea is 
described in the following proposition. 

Proposition 5.2. Let c = minx^dD 4>{x', p, p') and satisfy c > \p — p'\- Fix < s < rj'. 

If q E Aqd{p,p'), then Ago {p,p' + sA') = {q}. Moreover, if s < |(c— \p — p'\), then the 
operator Sg {Ec-s {p,p' + sA')) — Sg{dD) is positive definite on the common tangent space 
Tg{dD) = Tg{E,_sip,p' + sA')) = Tg{E,{p,p')) at q. 

Proof. Prom the definition of Aqd{p,p'), q G dD and one can write 

q^p' + s{A';p,p',c)A'. 

Then, we have 4>{q]p,p") = c — s, where p" = p' + sA'. Since (j){x;p,p") > c — s for all 
X G dD, this means c — s — minxedD (pi^'iP^p") and thus 

q e A9Dip,p' + sA'). 

Next let 

X G AQDip,p' + sA') . 
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We have (j){x]p,p") — minx^dD 4>{x; p, p") = c — s. Then we have 

c — s = 4>{x;p,p") = \p — x\ + \x — p"\ 

^\p-x\ + \{x-p') + {p' - p")\ 

> \p — x\ + \x — p'\ — \p' — p"\ = 4>{x;p,p') — s. 

Thus (f){x;p,p') < c and hence c = (l){x;p,p'). Prom these we have (l){x;p,p') — s — 
(f){x;p,p"). This is equivalent to 

\x — p"\ — \x — p'\ — s. (5.4) 

Since x is outside the open ball B" centered at p' with radius s, one can find the unique 
point x' on dB" such that jx — = miny^QB" \y — x\. Since we have |x— p'| = |a; — x'| +s, 
from (5.4) we obtain \x — p"\ = \x — x'\. Since p" e dB", it must hold that p" — x' and 
thus X = p' + s{A'; p, p', c)A' = q. 

The last statement is based on the fact that the eigenvectors for both shape operators 
are common and A' > A, where 

A = Xiq;p,p') = + A' = X{q;p,p' + sA') = + \ . 

\q ~ p\ \q — P\ IQ — P\ IQ—Pl—s 

See Appendix for these. Thus one concludes that the operator Sg {Ec-s {p,p' + sA')) — 
Sq{Ec{p,p')) is positive definite. Since Sq{Ec{p,p')) — Sq{dD) > 0, from these one gets 
the desired conclusion. Note that the condition s < ^{c—\p — p'\) is just for ensuring that 
\p - {p' + sA')\ <c-s. 

□ 

We have 

min (^(x; y, y') — min (f){x;p,p' + sA') — {rj + r]' — s) 

x£dD,yedB,y'edB^,_^{p'+sA') ^^^^ 

and, by (i) of Proposition 5.1, mmx^dD 4>ix; p, p' + sA') = c — s if g G Aqd{p,p')- 
Thus the condition 

T > min (l){x;y,y') 

xGdD,yGdB,y'edB^,_^{p'+sA ) 

is equivalent to (1.5). Therefore the following proposition is a direct consequence of 
Theorem 1.3 together with Proposition 5.2. 

Proposition 5.3. FixO < s < min {r]',{c-\p-p'\)/2). LetB and B' satisfy [BUB ]nD = 
0. Let f — Xb and g = Xn i a'\- ^ satisfy (1.5). Let q G Aqd{p,p') and set 

By/_g\P ) 

c = (t){q;p,p'). 

If D is admissible and dD is , then we have 

lim rS^^^-"-"') I / vjx- [ Wfdx] 

^-"^ \Jb Jb^,_^{p'+sA) J 

TT I 7] \ I 7] — S \ 1 



\q-p\J \\q-p'\-sj ^det{Sq{E,_,{p,p' + sA')) - Sq{dD)) ' 
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It is quite interesting to know the quantities contained in det {Sq{Ec-s{p,p' + sA')) — 
Sq{dD)). The following lemma whose proof is given in Appendix clarifies them. 

Lemma 5.1. Let < s < min(?7', (c — \p — p'|)/2). Let q e Aqd{p,p') and set c — 



det{Sg{E,_,{p,p' + sA')) - S,{dD)) = 



X{q;p,p' + sA') 



r\2 



jX{q;p,p' + sA') HQD{q) - , ^ + Kaoiq). 



1 + A- A 



2(1 + A- A') 



(5.6) 

Now we are ready to describe the proof of Theorem 1.4. 

Choose < si < S2 < min(?7', (c— \p — p'\) /2). Let s — si,S2. By (5.5) in Proposition 
5.3, one gets det {Sq{Ec-s{p,p')) — Sq{dD)) for s — Si, S2- Prom (5.6) we obtain the system 



/ 



2 \ 
jX{q;p,p' + SiA') 1 



1 + A • A' 



V V 1 + A-A 



jX{q;p,p' + S2A') 1 



(5.7) 



where 



V 



2(1 + A • A') 



I 



and 



^ det (5,(^C-..(P,P' + s,Al)) - Sq{dD)) - KQ;P,P' + siAr \ 



^ det {Sq{E,_,,{p,p' + S2A')) - Sq{dD)) 



X{q;p,p' + S2A) 



f\2 



I 



Since X{q]p,p' + SiA') < X{q;p,p' + S2A'), (5.7) is uniquely solvable with respect to X. 

This completes the proof of Theorem 1.4. 
Remeirk 5.2. It follows from (5.6) that 

det {Sg{E,{p,p')) - S,{dD)) - Paoil^; q) 



2 5q(9L>)(A X A') • (A X A') 



+ A • A' (1 + A • A') 



where 



^_ A(g;p,pO _ 1 1 ^ ^ 1 



2 2 \\q — p\ \q — p'l J 

Since this right-hand side has a bound 0{\A x A'p), this formula indicates an effect of 
the bistatic data on det {Sq{Ec{p,p')) — Sq{dD)). 
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6 Summary and some of open problems 



This paper is concerned with an inverse obstacle problem which employs the dynamical 
scattering data of acoustic wave over a finite time interval. The unknown obstacle D is 
assumed to be sound-soft one. The governing equation of the wave is given by the classical 
wave equation. The wave is generated by the initial data which is a characteristic function 
of an open ball B centered at p and observed over a finite time interval on a different ball 
B' centered at p' . It is assumed that [BVJB']r\D = 0. The observed data are the so-called 
bistatic data. This is a simple mathematical model of the data collection process using an 
acoustic wave/electromagnetic wave such as, bistatic active sonar, radar, etc. This paper 
aims at developing an enclosure method which employs the bistatic data. 

It is shown that from the data with some additional assumptions on the lower bound 
of T one can extract: 

(i) the first arrival time in the geometrical optics sense, that is, the shortest length of 
the broken paths connecting p to a point q e dD and q to p; 

(ii) the first reflection points between p and p', that is, all the points q e dD that 
minimize the length of the broken paths connecting p to q and q to p. 

(iii) the tangent planes of dD at all the first rcficction points. 

It is also shown that, under the admissibility condition for D, one can extract the 
Gauss curvature at an arbitrary first reflection point and the mean curvature with an 
additional term which depends on the positions of p, p' and the first reflection point. As 
a byproduct, for an example, a constructive proof of a uniqueness theorem for a spherical 
obstacle using the bistatic data is also given. 

We think that the problem taken up in this paper is a prototype of other various 
interesting problems. It is quite interesting whether the approach presented here can be 
applied to them or to develope its necessary modiflcation. Here we mention some of them. 

• Consider the sound-hard obstacle case or the obstacles with a dissipative boundary 
condition (cf. [21]). And also it is quite important to consider the corresponding problem 
for the Maxwell system. These remain open. 

• It would be interesting to consider also the case when obstacles are embedded in 
one of the two layers with known different propagation speeds and both the source and 
receivers arc placed in another layer. 

• Maybe the most interesting problem is that of extracting geometrical information 
about an unknown obstacle D behind a known obstacle Dq from the monostatic or bistatic 
data over a finite time interval. B and B' satisfy [BUB']n{DQUD) = and D is optically 
invisible from B and B' because of the existence of Dq. 

• It would be interesting to consider the case when B' is placed in the shadow region 
of D with respect to -B. It means that B' is invisible from B because of the existence of 
D between them. This is the case when [B U B'] fl dD ^ 0. What information about D 
can one extract from on i?'x ]0, T[? 

Although the author's interest is pursuit of the possibility of the enclosure method 
itself, research by other approaches to the problem taken up in this paper is also expected. 
Someone may think about the use of geometrical optics in the time domain as Majda 
has done in [19] for the problem considered in [18]. See also pages 440-447 in [25] for 
geometrical optics in the time domain and [22] for Majda's approach. His approach 
heavily depends on the hyperbolic nature of the governing equation in contrast to our 
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approach and it should be noted that the existing results by our approach can cover 
inverse problems for different type of equations like elliptic, parabolic or hyperbolic ones. 
Since this paper has not aimed at the comparative study of various approaches, we leave 
it to other opportunities. 
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7 Appendix 

7.1 Proof of (1.7) 

First we prove that 



mm (j){x;p,p') - {r] + 7]') > mm (j){x;y,y'). 

xeoD xedD,yedB,y'edB' 



Choose q e dD such that 



xGaD 



One can find yo G dB fl [q,p] and y^ G dB' n [q,p']. Since \yQ — q\ = \p — q\ — rj and 
\y'o ~ l\ — \p' ~ q\ ^ v'j we have 

yo, y'o) = p, p') + v') ■ 

This yields 

<l>{q;P,p') - {V + V')> min (f){x;y,y'). 

xGoD, yGoB, y'GoB' 

Next we prove that 

rmn(f){x;p,p')-{rj + r]') < mm (j){x;y,y'). 

xGaD xGaD,yGaB,y'£dB' 

Choose q e dD, yo e dB and y'g e dB' such that 

2/0, Z/o) = min d){x-y,y'). 

x&D, y&aB, y'eoB' 

Let y{s) be an arbitrary curve on dB such that y(0) = yo- Since (p{q',y{s),yQ) takes its 
minimum value at s = 0, we have 

^'P{q;y{s):yo)\s=o = o. 

This gives 

^■?(o) = o. 
\q - yo\ ds 
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Since dy/ds{0) can be an arbitrary vector perpendicular to the normal vector at yo e dB 
and q is outside of we have 



g - I/O ^ yo-p 
\(i-yo\ \yo-p\' 

This yields yo G [q,p] and similarly wc have y'g G [q,p']- Thus we have ||/o~g| = \p~q\ —V 
and ||/q — Q'I = |p' — ^1 — Tj'- This yields 

yo, y'o) = p> p') -{v + v') 

and thus 

0(?;yo,2/o) > mm - (rj + r]'). 

xGaD 

7.2 Proof of Proposition 1.1. 

Proof of (i). Since every point on dD is a point of support of we have (3.6). And 
from (3.6) we have e°(y^) > —Vf{y) for ddoiy) « 1 and y E D. Since t'/ > and C/ < 
(maximum principle), we have, for dgoiy) « 1 and y E D, |e°(y'')| < Vf{y). This yields 
the desired estimate. 

Let c = min^^zgn (f){x;p,p'). Assume that there exist two distinct points q and q' on 
^doiPiP')- Since D is convex, D becomes also convex. Thus every point on [q,q'] are in 
D. Choose an arbitrary point q" on [q, q'] \ {q, q'}. We have q" G D. Since both q and q' 
are on Ec{p,p') with c = miuj-gg^) (l){x;p,p'), it is clear that every point y on [q, q'] \ {q, q'} 
satisfies (f>{y',p,p') < c and thus (j){q";p,p') < c. Since ()){x;p,p') > c for all x G dD and 
g" G -D U it must hold that q" G -D. Since p in \ there exists a point g'" 
on [q",p] such that g'" G Then it is clear to have (f){q'"',p,p') < 4>{(l"]P,p') and thus 
4>{q"'',p,p') < c. However, since q'" G dD, we have 4>{q"'',p,p') > c. This is a contradiction. 
Proo/ of (a). Since is in the half space {x — q) ■ Uq < 0, it is easy to see that 
Sq{dD) < as the quadratic form on Tq{dD). On the other hand, wc have Sq{Ec{p,p')) 
is positive definite as the quadratic form on Tq{Ec{p,p')) = Tq{dD) (see Appendix). Thus 
Sq{Ec{p,p')) — Sq{dD) is positive definite as the quadratic form on the common tangent 
space and this yields (1.14). 

7.3 The shape operator for a spheroid 

Since c > |p — we have \p,p'] n Ec{p,p') = 0. Thus, for all x G Ec{p,p'), 

x-P ^ x-p' ^ ^ 
\x — p\ \x — p'\ 

and hence 

^ ^ x-P x-p' ^ ^ 
\x — p\ \x — p'\ 

Since 

V(j){x;p,p) = + -, 

\x — p\ \x—p'\ 
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we conclude that Ec{p,p') is a C°° surface and clearly compact. Let Ux denote the unit 
inward normal for x e Ecip,p'). We have 

V(/. 

= 



|V0| 



and note that 



These give the expression 



^2{l + A{x)-A'{x)) 



where 



A{x) = ^, A {x) 



Ix— p|' \x — p'\ 

Now we are ready to prove the following proposition. 

Proposition A.l. Let Sx denote the shape operator at x E Ec{p.p') with respect to 
which is the unit inward normal to Ec{p,p'). We have, for all v e T^EdpiP'), 

S.{v) = ^^"^ (h - Ia{x) ® A{x) - Ia'{x) ® A'{x)) V, {A.2) 

^2{1 + A{x)-A'{x)) V 2 2 J 

where 



\x — p\ \x — p'\ 

Proof. Let x — x{ai, (T2) be an equation for Ec{p,p') around q e Ec{p,p'). It means that 
x{0, 0) — q and 0(x((7i, a"2);p,p') = c. Since 

d . I X — p dx 



we have 



d ^ , . 1 dx X — p X — p dx 
-A(x) = 



daj \x — p\ daj \x — p\'^ \x — p\ daj 



\x 



1 Ox 
{I3 - A{x) (g) A{x)) 



\x — p\ daj 
This together with a corresponding expression for {d/daj)A'{x) gives 

-{A{x) + A'{x)) = X{x)Is - rA{x) (8) A{x) - -A'{x) ® A'{x) ^ 



daj \ \x — p\ \x — p'\ J daj 
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Moreover, since 



we have 

d 



{A{x) ■ A'{x)) 



\x — p\ daj \x — p'\ daj 



^(1 + A(.) . (a(.) . * ) - ^(1 + . A'i.)) (a'(.) . * ) 



we obtain 

d 



d 



1 9(7^- 



{A{x) ■ A'{x)) 



2(1 + A(x) • 
1 / 



2(1 + A(x) • A'(a;))V2 \ \x - p\ 
Thus one gets 

a A(x) + A'{x) A(x) + A!{x) (I dx , 1 dx ' 



^'^i ^1 + A{x) ■ A'{x) 2^1 + A{x) ■ A'{x) - P\ d(Jj \x - p'| da 

dx 



I ^ \Mx)h - I ^ , A{x) ® A{x) - -. — ^—-A{x) ® A!ix)\ 

^i + A{x)-A!{x) \ k-p'l J 



35 



and thus 



/ , , ., . — T77T 9 A(x) + A (x) 

9"} sjl + A{x) ■ A{x) 

1 1 (9t 1 1 i9r 

^ • - -A'(x)-^—^ ■ A'(x) 



|x — p| 2 |x — p'l 9(j 

--A{x)X{x)^ ■ A{x) - -A'{x)X{x)-^ ■ A'{x) + X{x) 



J 



2 daj 2 9(7j 9(7j 

Therefore we obtain 

d A{x) + A'{x) 



^1 + A(x) ■ A'(x) 

I ^=77== - ® M^) - \^{^) ® A\x)) ^ 

'l + A{x) ■ A{x) V 2 2 J daj 

and hence 

= , ^^'""^ , (h - Ia(x) ® A{x) - U'ix) ® 1^. 

^(^^^ ^2(1 + A(x) • V 2 2 ; daj 

Since 

c ( '^^ \ t \ \ — ( c ^ I I ^ 

\ dai c/(72 / \ oai da2 J 

we obtain (A. 2). 
□ 

Prom (A. 2) one can compute the principle curvatures at x e Ec{p,p'). 
Consider the case A{x) ^ A!{x). Since A{x^ and A'(a;) are unit vectors and A{x) ^ 
—A\x), we have v — A{x) x A!{x) ^ 0. Since A{x) ■ v — A!{x) ■ v — Q, v satisfies 

Sxyv) = I , V. 

^2{l + A{x)-A\x)) 
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Next choose v' — A{x) — A!{x). Since and A{x) + are parallel, v' e TxEc{p,p'). 

We have 

_ A(x)Vl + A(.).A-(a;) 

Note also that v and v' are perpendicular to each other. Therefore the eigenvalues of Sx 
consists of two real numbers: 



^2(1 + Aix)- A' (x))' 

iA.3) 

X{x)^l + A{x)-A'{x) 
'^^"^ ^ W=2 ■ 

If A{x) = A'{x), then i/^ = A{x). Since v • = for all v e Tx{Ecip,p'), from (A.2) 
we obtain 

Thus the set of all eigenvalues of consists of only \{x)/2. Therefore ki{x), k2{x) given 
by (A. 3) covers also this special case. Note that k2{x) < ki{x) and ki{x) — k2{x) if and 
only if A{x) = A'{x). Therefore the Gauss curvature K{x) at x e Ec(p,p') and the mean 
curvature H{x) with respect to are 

K{x) = ki{x)k2{x) = 

Hix) = h(^l±M^ - + ^(-) • ^'(-))- 

7.4 Proof of Lemma 5.1 

Replacing p' with p' + sA\ it follows from (4.18) and (4.19) that 



det {S,{E,^s{P:P' + sA')) - S,{dD)) = det ( M£lMl±i£)4 _ m,,- 

V ^2(1 + A.A') 

Xig;p,p' + sA')y A(,;p,p' + sA^^^^ ^^^^^ ^ ^^^^^^ 



2(1 + A ■ A') / J2{1 + A-A') 



From (4.14) we have 

ll + A-A' 



A-Ug^ A' -l/g^ 



{AA) 



where 

= Kq;p,p' + sA') ^^ . . ^ . . ) ^ ^(0). (A.5) 

2^2(1 + A- A') ^cTfcCT, 
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^(|A-efep+|A'-e,p) = l-A-A'. 

k=l 



TVace (M,,) = Kl:P,P' + ^A'}(1 -A- A') ^ ^^^^^^^^ ^^^^^ 



These give 



Thus one gets 



2^2(1 + A- A') 

For the computation of det (M^j) we prepare the following formula. 

Proposition A. 2. Let B be a 2 x 2-matrix, c and c' be two-dimensional vectors. Let 7 
be a constant. Let 

We have 



M = 7(c (8) c + c' (g) c') + B. 



(A.7) 



detM 



= 7^ det 



( ci c'l 



\ C2 C2 



+ 7 



B 



/ C2 



+ B 



( c' 



/ c' 



+ detB. 



Proof. We have 



detM 



(7(c? + cf ) + 611) (7(c^ + ) + ^22) - (7(ciC2 + c;c'2) + 6i2)(7(c2Ci + 4c;) + 621) 
= 7'(c? + cf )(c^ + + ^((c2 + cf )622 + (c^ + cf )6n) + 611&22 

-7^(ClC2 + c'iC2)(c2Ci + CaC'i) - 7((CiC2 + 44)621 + (C2C1 + 44)612) - 612621 

= f((c? + 4^)(c^ + 4^)-(ciC2 + 44f) 

+7((c^ + 4^)622 + {cl + 4^)611 - (ciC2 + 44)^21 - (C2C1 + 44)^12) + det B 

= 7^(ci4 - 4c2)^ 
+7((4 + 4')^'22 + (4 + 4')^'ii - (C1C2 + 44) (^'12 + &21)) + det B. 



□ 



Note that {Mtj) given by (A. 5) coincides with (A. 7) in the case when c — {A ■ e^), 
c' = (A'-ej), B^V^f{Q) and 



7 = 



X{q;p,p' + sA') 



2^/2(1 + A -A') 
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Prom (4.15) we have 

/ A • ei A' • ei \ 

det = det 

\ A - €2 A' ■ €2 / V A • 62 -A ■ 62 



A • 6i -A • 6i 

■ =0; 



from (4.19) we have 



vV(o) 



( A-e2 

-A ■ 6i 



A-62 \ 

-A-ei J 



= Sg{dD){{A ■ 62)61 - (A ■ 61)62) • ((A • 62)61 - (A • 61)62). 
By Lemma 4.3, we have 

A X A' = -^2(1 + A- A')((A • 62)61 - (A • 61)62). 

This gives 

S,{dD){{A ■ 62)61 - (A • 61)62) • ((A • 62)61 - (A • 61)62) 

1') ■ (A X A') 



one 



and thus from (A. 9) 

/A-62 

VV(0) 



Sq{dD){A X A^ 

2(1 + A- A'; 



V 



Since 



gets 

\ / A-62 
■A - 61 / V -A - 61 



Sg{dD){A X A') - (A X AQ 
2(1 + A - A') ■ 



we obtain also 



(A' - 62)61 - (A' - 61)62 = -((A - 62)61 - (A - 61)62), 



/ 


A' -62 \ 


/ 


A'-e2 \ 


V 


-A' - 61 J 


V 


-A' - 61 j 



S,(9i^)(A X A') • {A X A') 
-A- A') ■ 



l^-A'-6j V-A'-6j 2(1 + . 

Prom Proposition A.2, (A.8), (A. 10) and (A. 11) we obtain 

det (M,,) = -||^i^±^5,(ai))(A X A') - (A x A') + KMq)- 

Substituting this together with (A. 6) into (A. 4), we obtain (5.6). 
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